The Effect of the Environment on ct-Al 2 03 (0001) Surface Structures 
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We report that calculating the Gibbs free energy of the a-A^Os (0001) surfaces in equilibrium 
with a realistic environment containing both oxygen and hydrogen species is essential for obtaining 
theoretical predictions consistent with experimental observations. Using density-functional theory 
we find that even under conditions of high oxygen partial pressure, the metal terminated surface is 
surprisingly stable. An oxygen terminated Q-AI2O3 (0001) surface becomes stable only if hydrogen is 
present on the surface. In addition, including hydrogen on the surface resolves discrepancies between 
previous theoretical work and experimental results with respect to the magnitude and direction of 
surface relaxations. 

PACS numbers:68.35.Bs,68.35.Md,71.15.Ap 



The nature of the corundum surface (CV-AI2O3) is of 
considerable importance in a wide variety of technolog- 
ical applications. These range from catalytic supports 
and thin-film substrates to corrosion and wear protec- 
tion in mechanical systems. Yet, despite considerable 
experimental and theoretical efforts over the years [Q, 
the surface structure, and even the surface stoichiom- 
ctry, is a matter of strong controversy. The principal 
reason for this is that it is difficult to prepare clean, uni- 
form surfaces with specific, well-defined structures and 
stoichiometries. The typical heterogeneity of the surface 
makes the interpretation of low energy electron diffrac- 
tion (LEED) and other surface spectroscopic data diffi- 
cult. In addition, corundum is an insulator which makes 
the application of STM and other techniques based on 
electron spectroscopy problematic. 

Knowing the structure and stoichiometry of the corun- 
dum surface, however, is essential for understanding the 
electronic, mechanical, and chemical properties which de- 
termine its reactivity and performance in various applica- 
tions. Yet it has only been within the past two years that 
new experimental techniques have enabled the chemical 
identification of the surface terminations for these sys- 
tems. Renaud has reported a (1 x 1) structure, prepared 
by heating in an oxygen-rich atmosphere, which could 
reasonably be interpreted as being O-terminated [Q. 
Toofan and Watson found in a tensor LEED experiment 
both Al- and O-terminated domains in a 2:1 ratio, respec- 
tively ||. Ahn and Rabelais annealed the surface under 
UHV conditions and determined the detailed structure 
of the Al-terminated surface using time-of-flight scatter- 
ing and recoiling spectrometry (TOF-SARS) f|. The 
sensitivity of TOF-SARS also enabled the detection of 
hydrogen randomly distributed on the surface which was 
stable even at an annealing temperature of 1100° C. 

Despite recent advances in experimental techniques, 
however, many questions and inconsistencies remain that 
theoretical calculations performed to date have not been 



able to resolve. Corundum (a-A^Os) crystallizes in a 
structure which can be described by a primitive rhombo- 
hedral unit cell with two AI2O3 formula units (10 atoms) 
or by a conventional hexagonal unit cell with six AI2O3 
formula units (30 atoms). For the hexagonal unit cell, the 
atoms are stacked along the (0001) direction according to 
the sequence i2-AlA103-AlA103-i2, where R represents 
the continuing sequence in the bulk. The oxygen atoms 
form hep layers, and the metal atoms fill two thirds of 
the octahedral sites between these layers || . 

In the i?-Al-Al-03-i? corundum structure (0001) stack- 
ing sequence, there are three unique stoichiometric slice 
planes. Yet theoretical methods ranging from the empir- 
ical to the ab initio P [Tl|] have so far identified only one 
stable termination stoichiometry for the a-A^O^ (0001) 
surface: an Al-monolayer indicated by AlOsAl-i? in Fig. 
1. No stable O-terminated surface has been found, so 
how can the experimentally observed O-terminations be 
explained? In addition to the surface stoichiometry, 
there is considerable disagreement with respect to the 
magnitude and direction of the surface relaxations. For 
the Al-terminated surface, experimentally observed re- 
laxations range from -51% @ to -63% (§, whereas the 
density-functional theory (DFT) pseudopotential calcu- 
lations predict a relaxation of about -87% ||]jn]], and 
Hartree-Fock gives -40% 0. Toofan and Watson re- 
ported an outward relaxation for their two domain sys- 
tem, which is the opposite of what other experimental 
and theoretical investigations had concluded. These dis- 
crepancies are significant and need to be resolved. 

What has not been previously considered in the the- 
oretical treatment of the CV-AI2O3 (0001) surface is the 
effect of the environment on the surface structures and 
stoichiometry. Under realistic conditions, a surface will 
exchange atoms with its surroundings. Hence in this 
work we present the first analysis of the Gibbs free en- 
ergy of the system with respect to its dependence on the 
chemical potentials of the components present in the ma- 
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terial and the environment at OK and 1000K. For a metal 
oxide such as a-A^Os, the O2 partial pressure is obvi- 
ously the most important factor in the analysis, as well as 
temperature. In addition, the presence of stable hydro- 
gen on the surface also needs to be addressed, as it can 
be incorporated into the bulk structure during growth, 
remain from calcination of Al(OH)3 in the synthesis of 
a-Al203, or may result from exposure of the surface to 
water vapor prior to placement in the UHV chamber. 
Yet to our knowledge the presence of hydrogen on the 
corundum surface has not yet been investigated in theo- 
retical studies. In this paper we report that the surface 
stoichiometries, structures, and properties change signif- 
icantly depending upon the chemical potential of O2, H2, 
and H2O gases in equilibrium with the surface at differ- 
ent temperatures. We use the full-potential linearized 
augmented planewave (FP-LAPW) method to solve the 
Kohn-Sham equations and calculate the total energies, 
forces, and chemical potentials for all reasonable (1 x 1) 
corundum (0001) surface geometries. We find that the 
metal-terminated surface is surprisingly stable across the 
range of a physically realistic oxygen chemical potential. 
An oxygen terminated surface becomes stable only if hy- 
drogen is present on the surface, even at partial O2 pres- 
sures which are too high for standard UHV equipment. 
In addition, including hydrogen on the surface of both 
the aluminum and oxygen terminated surfaces results in 
calculated relaxations which agree with the latest exper- 
imental results, in both magnitude and direction. 

To systematically investigate the (lxl) (0001) sur- 
face which is observed at annealing temperatures below 
1250° C, we have generated what we believe are all the 
possible (1 x 1) (0001) Al- and/or O-surface terminations 
of AI2O3, plus several additional surfaces containing hy- 
droxyl groups. Previous studies which investigated possi- 
ble (1 x 1) (0001) surface terminations examined only the 
three types of structures obtained by simply cleaving the 
unit cell at unique positions [p|JIo|: A10 3 Al-i?, 3 AlAl-i?, 
and AIAIO3--R. For the single Al-terminated layer there 
are three other possible locations for the Al atom in ad- 
dition to the corresponding bulk site, as shown in Fig. 
For the O-terminated surface, other investigations have 
not considered the possibility of oxygen vacancies which 
could occur and still maintain a perfect (lxl) surface 
periodicity. Two additional O-terminated structures can 
be created by introducing one or two oxygen vacancies 
per unit cell, indicated by C^AlAl-i? and OiAlAl-i?, re- 
spectively. Hence, there are a total of at least eight pos- 
sible (1 x 1) geometries of the (0001) surface which we 
investigate to determine their relative stability. 

In our calculations, the (1 x 1) CV-AI2O3 surface is mod- 
eled by a slab which consists of a finite number of layers 
and is infinite in the plane of the surface. The slabs are 
repeated periodically along the [0001] direction and sep- 
arated by 10A of vacuum. The slab contains six oxygen 
O3 layers and from ten to fourteen aluminum layers, de- 




FIG. 1. Side and top view of the Al-03-Al-i? surface of 
a-Al2C>3 , in which 1, 2 and 3 indicate three sites which could 
be occupied by the topmost Al instead of the bulk site. The 
large circles denote O atoms and the small circles Al atoms. 



pending upon the specific surface studied. We carefully 
tested that the thickness of the vacuum as well as that 
of the slabs are sufficiently large to ensure that surface- 
surface interactions through both the vacuum and the 
slab are negligible. jl3| The two surfaces of the slab are 
identical and inversion symmetry is maintained. 

For our total-energy calculations we use the general- 
ized gradient approximation of Perdew et al. ]T^ | for 
the exchange-correlation potential and the FP-LAPW 
method as implemented within the WIEN97 program 
Jl5| [TtJ to solve the Kohn-Sham equations A uni- 

form k-point mesh with ten points is used for the entire 
surface Brillouin zone. An identical mesh is used for the 
bulk calculations to ensure consistency. The calculated 
equilibrium lattice constants for the hexagonal unit cell 
and the non-symmetry fixed positions of the atoms agree 
with experimental values to within 0.7%. The calcula- 
tions give the heat of formation at OK of the bulk AI2O3, 
AH°— 17.37 eV, which is in good agreement with the ex- 
perimental value of 17.24 eV pl|. 

To compare the stability of surfaces with different sto- 
ichiometries in chemical and thermal equilibrium with 
the gas phase and with the bulk, we calculate the Gibbs 
free energy of the surface relative to the chemical po- 
tential of oxygen. The Gibbs free surface energy of a 
slab at temperature T and partial pressure p is given by 
ft = S total + AG vib -X;^^(T,p), where £ total is the scf 
energy of the slab, AG vlb is the vibrational contribution 
to the Gibbs free energy, Ni the number of zth type of 
atoms in the slab, and fii(T,p) is the chemical potential 
of the ith type of atom at a given temperature and pres- 
sure. The sum is over all the types of the atoms in the 
slab. To calculate the values of the chemical potential for 
all species we use /j,i(T,p) = $ + A^(T,p), where OK is 
taken as the reference state, and A/Xj(T,p) is the change 
in free energy from that reference state to the system 
at a given temperature and pressure. The OK values for 
surfaces and for bulk Al and AI2O3 are taken from our 
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total energy calculations. For the dissociation energies 
of the H2, O2, and H2O molecules, we use experimental 
values. The A/Zj(T,p) values for all crystalline and gas 
phase species are taken from the JANAF Thermochemi- 
cal Tables |o). 

As each surface is considered here to be in chemical 
and thermal equilibrium with the bulk and the environ- 
ment, the chemical potentials are constrained, so for the 
clean surface: 2/zai + 3^o = ^Ab0 3 ' Mai < E^ lk , and 

Mo < ^o 2 ° lcCulC : where e ai'o 3 ^ ft** 3 total energy per 
bulk AI2O3 formula unit. The limits of the chemical po- 
tential for oxygen are determined by conditions of equilib- 
rium with relevant oxygen-containing species in the sys- 
tem. The maximum oxygen chemical potential is that of 
a maximum concentration of O2 molecules at OK, which 
corresponds to O2 condensing on the surface. This is the 
reference value of the oxygen chemical potential, which 
is the zero on the right hand side of the graph in Fig. 
^. Since the system is in equilibrium with bulk AI2O3 
the minimum no occurs when /xai is at a maximum and 
fio = ^M% 3 ~ ^M lk - Below this range, aluminum 
metal would condense on the surface. 

For the surfaces where hydrogen is adsorbed, the free 
energy of the surface is calculated relative to the chemical 
potential of hydrogen /in in both H2 and H2O, as desorp- 
tion of hydrogen as either molecular species is possible. 
The lines for the H0 3 AlAl-i? and HsOsAlAl-i? surface 
energies calculated with respect to fin in equilibrium with 
H 2 and H 2 cross at no equal to -2.7 eV at OK. Below 
-2.7 eV, /ih is determined by the equilibrium with H2, 
and above with H2O. Below these lines either H2 or H2O 
would condense on the surface. It should be noted that 
the only independent variables in our system is no m 
equilibrium with bulk AI2O3, and nil- From these the 
dependent Al, O2, H2 and H2O chemical potentials fol- 
low. 

The results are shown in Fig. |2|. The Gibbs free en- 
ergies per surface area for various (lxl) geometries of 
the (0001) surface are displayed as a function of no- The 
dashed black vertical lines bracket the allowed range of 
no- For the systems containing just aluminum and oxy- 
gen, there is clearly only one overwhelmingly low energy 
surface stoichiometry: the AlO^Al-i? structure. This is 
consistent with the previous theoretical results. Our re- 
sults show, however, that this stability extends across the 
entire range of the oxygen chemical potential, even up to 
the limit where an O2 condensate would form on the 
surface. The empirical models explain this remarkable 
stability on the basis of simple ionic considerations, as it 
is the only structure without a significant surface dipole. 
The actual situation is a bit more complicated, of course, 
and the quantum mechanical calculations performed by 
us and others I^JlCfl indicate that the surface aluminum 
atom relaxes inward until it is practically coplanar with 
the oxygen layer, and rehybridizes to an sp 2 orbital con- 
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FIG. 2. Surface energies of different AhC^OOOl) surface 
terminations. The dashed verticle lines indicate the allowed 
range of the oxygen chemical potential, no- Values in red 
indicate examples of higher temperatures and pressures. The 
green and blue regions indicate the range where hydrogen on 
the surface is in equilibrium with H2 and H2O, respectively, 
from OK up to 1000K and 1 atm pressure (red lines). 

figuration which is dramatically stabilized by autocom- 
pensation and bonding considerations. That is, the struc- 
ture relaxes and charge is transferred so that the sp 2 
bonds between the aluminum atoms and the three oxy- 
gens are filled, the aluminum 3p z orbital perpendicular 
to the surface is empty, and there are no partial occupan- 
cies of dangling bonds at the surface. In our calculations 
the lowest empty surface state, which is primarily the Al 
3p z with some Al 3s and O 2p character, is 4.5 eV above 
the Fermi level. This is very high and contributes to the 
relative chemical inertness of the relaxed Al-terminated 
surface. Moving the topmost aluminum atom of the Al- 
tcrminated surface from its position to one of the three 
non-bulk positions indicated in Fig. [l], increases the sur- 
face energy dramatically from 133 meV/A 2 to 355, 282, 
and 228 meV/A 2 for the three sites, respectively. 

The 03-terminated surface has a very unfavorable en- 
ergy, although this decreases sharply as no increases. 
This stoichiometry has the highest dipole at the surface, 
which is reflected in a very large work function of 9.66 eV, 
compared to 5.97 eV for the Al-terminated surface above. 
In addition, there are partially filled open valencies on 
the surface oxygens which are not relieved by relaxation. 
Our results indicate, however, that this high energy sur- 
face can become more stable through oxygen evapora- 
tion. The 02-terminated surface has a consistently lower 
free energy than the O3AIAI-J?, and the OiAlAl-i? sur- 
face is more stable than O2A1AI-.R, except at the highest 
range of the oxygen chemical potential. At 1000K the 
metal-terminated surfaces become the lowest energy sur- 
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TABLE I. Inter layer relaxations at the Al-terminated sur- 
face in % of the corresponding bulk spacings. 



interlayer 


Theory of Ref. 
this § @ 


Experiment of Ref. 
1 0] 1 


A1-O3 1-2 


-86 


-86 


-87 


-63 -51 +30 


O3-AI 2-3 


+6 


+3 


+3 


+16 +6 


Al-Al 3-4 


-49 


-54 


-42 


-29 -55 


AI-O3 4-5 


+22 


+25 


+19 


+20 


O3-AI 5-6 


+6 




+6 





face under the most oxygen-deficient conditions. This is 
consistent with what happens at annealing temperatures 
above 1250° C, where oxygen evaporation leaves behind 
a metallic aluminum overlayer [jj]]. The figure does not 
convey the kinetic barriers required to move from one 
surface stoichiometry to another, but the relative ther- 
modynamic stability is clear and is consistent with ex- 
perimental observations. 

The large surface dipole and dangling bonds of the 
03-terminated surface can also be compensated by the 
addition of hydrogen to the surface. Adding one hydro- 
gen per unit cell to the surface |2lJ dramatically low- 
ers the free energy and the work function from 9.66 to 
7.06 eV. This surface is much more stable than the 03- 
terminated surface at all oxygen partial pressures. Sat- 
urating all surface oxygens with hydrogen, indicated by 
the HaOsAlAl-i? line in Fig. 2, results in the lowest free 
energy and greatest stability across the entire range of 
physically realistic conditions. The oxygen open valen- 
cies are filled, charge transfer occurs from the H to the 
topmost O layer, and the work function is reduced to 
3.44 eV, the lowest of all surfaces examined. The region 
of negative surface energy of the H3 03-terminated sur- 
face indicates that the hydroxylated surface is lower in 
energy than bulk ev-Al203 plus free H2O. This negative 
surface energy reflects the strength of the H-OA1 bond 
and is consistent with the negative heat of the CV-AI2O3 
+ 3H2O — > 2A1(0H)3 reaction. This reaction is exother- 
mic by 0.19 eV per AI2O3 formula unit relative to liquid 
water at 298. 15K and 1 atm pressure f22|| . 

We find that in addition to explaining the relative sta- 
bilities of the (0001) (1 x 1) surfaces observed experi- 
mentally, considering the presence of hydrogen on the 
surface is essential to resolve the controversy surrounding 
the surface relaxations. Our results are shown in Table 
H In the absence of hydrogen, our full potential results 
find a large inward relaxation of 86% for the first layer, 
in very close agreement with previous pseudopotential 
calculations although there is less agreement for 

subsequent layers. When a hydrogen atom occupies site 
1, however, the contraction between the first Al and O 
layers is ~ 69%, in very close agreement with the ex- 
perimental result of Ahn and Rabalais Q . Several other 
locations for the hydrogen atom were investigated, with 
similar results. 



For the oxygen-only terminated surface, theoretical 
predictions have indicated an inward relaxation of the 
top layer of 0.05A. Yet the experimental results of Toofan 
and Watson, in which one of the two domains observed 
was O-terminated, clearly showed an outward relaxation 
of 0.12A ||. Placing a hydrogen at site 1, however, re- 
sults in a predicted outward relaxation of 0.11 A, in excel- 
lent agreement with experiment. Saturating the surface 
with hydrogen, the H 3 03-structure, increases the pre- 
dicted relaxation to 0.19A. Hence even though H was not 
detectable with the instrumentation used by Toofan and 
Watson, our free energy calculations and predicted re- 
laxations provide strong evidence that H was most likely 
present on the surface. 

Real surfaces are usually much more complicated and 
heterogeneous than the model (lxl) system studied here, 
but encorporating the effects of chemical and thermal 
equilibrium with the environment into the model pro- 
vides a sound framework for interpretation of experimen- 
tal results on these real systems. 
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